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Abstract 

Asymmetric nuclear matter equation of state at finite temperature is studied in 
SU(2) chiral sigma model using mean field approximation. The effect of temperature 
on effective mass, entropy, and binding energy is discussed. Treating the system as one 
with two conserved charges the liquid-gas phase transition is investigated. We have 
also discussed the effect of proton fraction on critical temperature with and without 
p-meson contribution. We have extended our work to study the structure of proto- 
neutron star with neutron free charge-neutral matter in beta-equilibrium. We found 
that the mass and radius of the star decreases as it cools from the entropy per baryon 
S = 2 to S = and the maximum temperature of the core of the star is about 62 MeV 
for S = 2. 
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1 Introduction 



The study of properties of hot dense asymmetric nuclear matter has been, in the last 
few years, vigorously pursued in connection with astrophysical problems[l,2], such as 
supernova explosions, the evolution of neutron stars and nuclear problems such as 
heavy-ion-coUisions. As the equation of state(EOS) describes the variation of energy 
density and pressure with density and temperature, it can be used to describe different 
phases like gaseous and liquid nuclear matter upto the deconfinement transition. It is 
also possible to study the liquid-gas phase transition, which may occur in the warm 
and dilute nuclear matter produced in heavy-ion collisions. Several authors using 
non-relativistic[3] and relativistic[4-8] theories have studied hquid-gas phase transition. 
Most of the calculations found the critical temperature Tc lying in the range of 14-20 
MeV for symmetric nuclear matter. The critical temperature of symmetric nuclear 
matter for the most acceptable relativistic Walecka model, is Tc ~ 18.3 MeV[9]. As the 
asymmetric parameter or the proton fraction plays a vital role in getting the critical 
temperature, the addition of p-meson is quite essential for the study of asymmetric 
nuclear matter [4, 5]. 

Field theoretical finite temperature EOS plays an important role in studying the 
properties and structure of hot dense massive stars such as proto-neutron stars at dif- 
ferent temperatures and densities. The important characteristics which determine the 
composition of matter in a compact star are[10] their relative compressibilities (important 
to determine maximum mass of neutron star), symmetry energies ( important to de- 
termine the typical stellar radius and the relative n, p, e, neutrino abundances) and 
specific heats (important to determine local temperatures). These characteristics vary 
with respect to the EOS used in different models. Generally the structures of both 
hot and cold, and both neutron- rich and neutron-poor, stars are fixed by the EOS [2]. 
A Proto-Neutron Star(PNS) is born following the gravitational collapse of the core of 
a massive star during a supernova explosion (type-II) and evolves to a cold and delep- 
tonized neutron star, basically through neutrino emission. This very dense and hot 
core is also able to trap neutrinos, imparting momentum to the outer layers and then 



cooling as it reaches a quasi-equilibrium state. There can also be a quark-hadron phase 
transition in PNS at high density and temperature[ll]. 

In this paper we have extended our earlier study[12] on finite temperature EOS 
by adding p-meson to the Lagrangian density within the MCH model [13] as the ad- 
dition of p-mcson is quite essential for the study of asymmetric nuclear matter. In 
our earlier work we have studied the effect of temperature on EOS, effective mass, en- 
tropy and binding energy for symmetric nuclear matter and investigated the liquid-gas 
phase transition using the same model without the p-meson. We now study here how 
the asymmetric parameter or proton fraction in addition to p-meson changes those 
properties at finite temperatures. As finite temperature EOS has an important role 
in studying the structure and properties of astrophysical objects, we have also inves- 
tigated the structure of PNS taking /^-stable charge neutral matter without neutrinos 
consisting of neutrons, protons and electrons only. 

In Sec. 2, asymmetric nuclear matter at finite temperature is presented. Proto- 
neutron star is discussed in Sec. 3. We conclude with some remarks in Sec. 4. 



2 Asymmetric nuclear matter at finite temperature 
2.1 Equation of state 

Continuing our earlier investigation [12] on the study of finite temperature effect of 
asymmetric nuclear matter using MCH model we include here the p-field in the La- 
grangian density. The EOS for hadronic phase is calculated by using the Lagrangian 

density [13] (with h = C= KBoltzmann = 1), 

AC 

OTTl 

In the above Lagrangian, F^^ = d^u^, - d^oJ^,, G^^ = d^p^ - d^p^, and x = (tt^ + a^)^/^, 
-0 is the nucleon isospin doublet, tt is the pseudoscalar-isovector pion field, a is the 



scalar field and a;^, is a dynamically generated isoscalar vector field, which couples 
to the conserved baryonic current = '07/i'0- is the isotriplet vector meson field 
with mass nip. B and C are constant coefficients associated with the higher order 
self-interactions of the scalar field . 

The masses of the nucleon, the scalar meson and the vector meson are respectively 
given by 

m = g^XQ, rria = V2\xq, = q^xq (2) 

Here xq is the vacuum expectation value of the a field , g^j, Qp and are the coupling 
constants for the vector and scalar fields respectively and A = (m^ — ml)/{2f^), where 
is the pion mass , f^^ is the pion decay coupling constant . 

Using Mean-field approximation, the equation of motion for isoscalar vector meson 
field is 

riB 

uJo = ^ (3 

and that for p-field is given by 

pI = {gp/2ml){np - Tin) (4) 
nsi— np + Un) is the baryon number density at temperature T and is given by 

riB = ttP^ d^k [miT) - fiiiT)] (5) 



(27r)= 



with 



where i = n,p. E*{k) = (A;^ + y'^m'^y^'^ is the effective nucleon energy, P = 1/KbT, 
7 is the spin degeneracy factor with nj(T) and ni{T) being Fermi-Dirac distribution 
functions for particle and anti-particle respectively at finite temperature. Ui is the 
effective baryon chemical potential which is related to the chemical potential /li as 

^P^I^P —2 jK^p-^n) (8) 

Ct^UB Cp 

l^n=IJ'n — + ^{'np-nn) (9) 



and y = x/xq is the effective mass factor whicli must be determined self consistently 
from the equation of motion for scalar field which is given by 



0. 



(10) 

m* = ym is the effective mass of the nucleon and the coupling constants are 

2 2 

a = ^, = % and C, = gl/ml . (11) 

The symmetric energy coefficient that follows from the semi-empirical nuclear mass 
formula isflSl dsvm = i^-r H ^4^ ? ■ 

L J sym i27r2 6^kj+m*^ 

Now the nucleon number densities, energy density and pressure at finite tempera- 
ture and finite density are given by 

^ (2^ r "^'^ f""^^^^ " ^^^^^^ ^^^^ 

rin = ^ 0?'^ KiT) - nn{T)] (13) 

+ ^ c^fcA;V(A;2 + m*2) [n^(T) + ri„(r) + 71^(7) + np{T)\ (14) 
^ _m^(l-y^)^ _B_ C _ Ozl 

The entropy density(S/V) and entropy per baryon(S) can be obtained as 

S/V = (P e - jJipUp - iJinnn)P (16) 

S ^ {P + e - i^pUp - iJ,nnn)f3/nB (17) 
In order to describe the asymmetric nuclear matter one can introduce the proton frac- 



tion which is defined as 



yp-^ (18) 



For the neutron matter yp = and for symmetric nuclear matter i/p = 0.5. 

The values of four parameters C^, C^, Cp, B and C occurring in the above equations 
are obtained[13] by fitting with the saturation values of binding energy/nucleon(-16.3 
MeV), the saturation density (0.153 fm~^), the symmetric energy (32 MeV), the effec- 
tive(Landau) mass(0.85 M), and nuclear incomprcssibility(~300 McV), arc C^^ = 1.999 
fm2, Ca = 6.816 fm^, Cp = 4.661 fm^, B = -99.985 and C = -132.246. For a given value 
of ub at fixed i/p and/or T, the equations[10,12,13] can be solved self-consistently to 
get y, /ip, /In and making substitution, the energy density, pressure, entropy density 
and entropy per baryon can be evaluated. 

Now let us discuss the liquid-gas phase transition. In case of the more common 
single component phase equilibria, such as liquid vapour, the phases are distinguished 
by only one parameter, e.g. the density whereas the binary mixture has the additional 
parameter, the proton fraction y^, which is different from one phase to the other [14]. For 
the description of liquid-gas phase transition we have followed here the thermodynamic 
approach of Refs. [7,14,15]. In the case of asymmetric nuclear matter the system is 
characterized by two conserved charges namely the baryon density (n^ = Hp + n„) 
and the total charge or equivalently the third component of isospin(73 = "'^~"'" ) or 
equivalently the proton fraction (yp). Thus the stability criteria may be expressed by 
the following relations [7] 



The first inequality shows that isothermal compressibility is positive which implies 
that the system is mechanically stable. The second condition expresses the diffusive 
stabihty. If any of the three stability criteria is violated then there will be a phase 
separation [7]. 




(19) 




(20) 



2.2 Results and discussions 

First, we discuss the effective nucleon mass of nuclear matter at finite temperature and 
density. Fig.l sliows tlie effective mass versus number density at temperatures T=0, 
100, 200 MeV. Tlie solid lines are for yp = 0.5 and dotted lines are for yp = 0. It 
is observed that in all temperatures M* decreases with the increase of n^. For T=0 
or T=100 MeV, the effective mass varies little with yp, whereas it is very sensitive 
with yp for T:=200 MeV. It is also clear that at zero density, the effective mass is the 
same for neutron matter and symmetric matter at different temperature. Hence these 
result indicate that the nucleons of neutron matter must stay in higher energy levels 
compared to that of symmetric matter in order to have the same number density and 
the difference will decrease as the number density decreases. This result agrees with 
that obtained in Ref.6. 

In Fig. 2 we show the effective mass as a function temperature at zero density(nB = 
0) and saturation density (n^ — 0.153fm~^). The solid lines are for symmetric nuclear 
matter(yp=0.5) and doted lines for neutron matter(yp = 0). It is clear from figure that 
for ub — 0, the solid line and the dotted line are same, whereas for ub — tiq, the yp 
dependence is sensitive to the temperature in between 150 MeV and 240 MeV and again 
the two lines coincide at higher temperatures. It is also observed that for ub = no, 
the M* first increases slowly and then falls suddenly at about T fti 240 MeV. But for 
Ub — 0, M* remains almost constant as the temperature increases and falls suddenly 
at about T 235 MeV. This result shows that a first order phase transition appears at 
ns = 0, T Ri 235 MeV which is similar to the result obtained for Walecka model, which 
has such a phase transition at T ^ 185 MeV[16]. Because of strong attraction between 
the nucleons at high temperatures the nucleon antinucleon pairs can be formed which 
may lead to abrupt change in M* to take place in high temperature region. But the 
mechanism of this first order phase transition is not clear. It is expected that this 
abrupt change of M* in the high temperature region for some models including the 
Walecka model might be related to the formation of new matter [17]. 

Fig. 3 shows the binding energy per nucleon as a function of the baryon density at 
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Figure 1: Effective mass as function of baryon number density at different temperatures. 

different temperatures for symmetric nuclear matter, yp=0.5. At zero temperature it 
lias a minimum at tlie nuclear saturation density tiq which corresponds to a binding 
energy per nucleon of -16.3 MeV. With the increase of temperature the minimum 
shifts towards higher densities and for higher temperatures the minimum of the curve 
becomes positive. It is also observed that as the temperature increases, the nuclear 
matter becomes less bound and the saturation curves in the MCH model are seems to 
be flatter than those observed in Walecka model[17,18]. This result implies that the 
nuclear matter EOS in MCH model is softer than that obtained in Walecka model. 

In Fig. 4, we show pressure as a function of density ub at different temperatures. 
The solid lines are for symmetric nuclear matter and the dotted are for pure neutron 
matter which are stiffer at all temperatures. For a given tib, the pressure has the 
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Figure 2: Effective mass as function of temperature for constant baryon number densi- 
ties. 

usual trend of increasing with temperature [19]. As the temperature increases the EOS 
becomes stiffer. Pressure has a non-zero value for = at and above a temperature 
of 200 MeV. It indicates that pressure has a contribution arising from the thermal 
distribution functions for baryons and anti-baryons as well as from the non-zero value 
of the scalar field. Similar results were also obtained by Panda et al. [20] for symmetric 
nuclear matter. The non-zero value for scalar a field has also been observed in Walecka 
model[19]. 

The entropy density as a function of density at different temperatures for symmetric 
nuclear matter and pure neutron matter is presented in Fig. 5. It is observed that 
entropy density for both is non-zero even at vanishing baryon density at a temperature 
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Figure 3: Binding energy per nucleon as function of baryon number density at different 
temperatures. 

of 200 MeV with contributions from the non-zero value of the sigma field. Similar 
behavior was also observed for entropy density in the Walecka model and ZM model 
calculations [18] . This increase of entropy density with increase of temperature indicates 
a phase transition. The Entropy per baryon(S) as a function of baryon number density 
at different temperatures for symmetric nuclear matter is shown in Fig. 6. At lower 
temperatures S decreases slowly as compared to higher temperatures and the minimum 
value of S increases as the temperature increases which is similar to the result obtained 
in Refs.[2,8]. 

We now discuss the liquid- gas phase transition. The pressure as a function of baryon 
density at fixed temperature T = 10 MeV with different proton fractions is shown in 
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Figure 4: Pressure (P) as function of baryon number density at different temperatures. 

Fig. 7. It may be observed from the figure that for any fixed density with fixed T = 10 
MeV, the pressure is not constant rather it increases with decrease of proton fraction. 
This clearly indicates that for asymmetric nuclear matter during isotherm liquid-gas 
phase transition the pressure can not remain constant but increases monotonically. It 
shows that for small yp, particularly for neutron matter (^^=0), the pressure increases 
monotonically which indicates that matter is stable at all densities. But for yp > 0.2, 
the compressibility becomes negative, indicating mechanical instability. The diffusive 
unstable regions can be seen clearly from Fig. 8, where the chemical potentials of proton 
and neutron is shown as a function of yp at fixed pressure P = 0.1 MeV/fm^ and 
temperature T = 10 McV. According to the inequality[20] the region of negative slope 
for Hp and positive slope for /x„ is unstable. Thus violation of stability criteria is an 
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Figure 5: Entropy density as function of baryon number density at different tempera- 
tures. 

indication of phase separation. 

Fig. 9 shows the variation of pressure as a function of baryon density for different 
yp. One can see that the region of mechanical instabihty decreases both with increase 
of temperature and decrease of proton fraction [8]. We present the pressure of sym- 
metric nuclear matter as a function of baryon density for low temperatures in Fig. 10. 
The figure shows that at zero temperature, the pressure first decreases, then increases 
and passes through P = at = no (saturation density), where the binding en- 
ergy per nucleon is a minimum. Decrease of pressure with density implies a negative 
incompressibility, K = 9(^^), which is a sign of mechanical instability. When the 
temperature increases the region of mechanical instability decreases and disappears 
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Figure 6: Entropy per baryon as a function of baryon number density at different 
temperatures. 

at the critical temperature Tc, which is determined by 1^^= \tc— 0, above 
which the hquid-gas phase transition is continuous. We have obtained the value of 
critical temperature Tc ^ 17.2MeV, critical density ric ~ 0.045/m~^, critical pres- 
sure Pc ~ 0.27AMeV/ fm~^ and critical effective mass M* ^ 887 MeV which is in fair 
agreement with the results obtained in other studies [8, 17, 18]. 

In Fig. 11, we plot the variation of critical temperature with different proton fraction(yp) 
with and without p. The critical temperature Tc decreases monotonically[8,21] as the 
proton fraction decreases and goes to zero for yp — 0.02 with p whereas Tc — 11.6 at 
yp = without p. This indicates that addition of p-mcson lowers the critical tempera- 
ture at smaller yp. As may be seen from equations[8,9], the addition of p- meson gives a 
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Figure 7: Pressure as a function of baryon number density for different proton fractions 
at temperature T—10 MeV. 

repulsive potential which makes the nuclear matter easier to be gasified. But for neu- 
tron matter(yp = 0), the system only remain in gas phase even at zero temperature[6] . 

3 Proto-Neutron Star 

In this section we extend our investigation to a study of the structure and properties of 
Proto-Neutron Star. A PNS is born following the gravitational collapse of the core of 
a massive star during a Supernova explosion(type-II). It is a hot collapsed core which 
can reach temperatures as high as few tens of MeV. The evolution of PNS proceeding 
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Figure 8: Chemical potentials as a function of Up at temperature T=10 MeV and P=0.1 
MeV/fm^ 

through several distinct states with various outcomes is discussed in Ref.[2]. During 
the early evolution of PNS, a neutron star with an entropy per baryon of order of unity 
contains neutrinos that are trapped in matter on dynamical time scale and after a lapse 
of few tens of seconds the star achieves it's cold catalysed structure with essentially 
zero temperature and zero trapped neutrinos. A PNS has approximately uniform 
entropy per baryon(S) of 1-2 across the star [22]. At birth the PNS has S=l. After 
deleptonization the entropy per baryon reaches its maximum (S~ 2) and finally cools 
down to its cold state with S=0[2]. The finite temperature aspect of EOS plays an 
important role in the study of properties and structure of PNS. 

The structure of PNS mainly depends on it's composition [2]. Since the composition 
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Figure 9: Pressure as a function of baryon number density for different proton fractions 

of neutron star basically depends on the nature of strong interactions which are not well 
understood in dense matter, one has to investigate various possible conditions taking 
many possible models. Out of various possible cases discussed in Ref.[2], we consider 
here a case in which matter consists of neutrons, protons and electrons whose relative 
concentrations is determined from the conditions of charge neutrality and /3-equilibrium 
in the absence of neutrino trapping [2]. 

The /9-equilibrium (without neutrino trapping) and the charge neutrality conditions 
are respectively given by 

/^n = /Wp + fJ-e (21) 

and 

rip = He (22) 
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Figure 10: Pressure (P) as function of baryon number density for symmetric matter at 
different temperatures. 

Where /le and rig are chemical potentials and number density of electron respectively. 
The electron number density at finite temperature can be written as 



Tie. 
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d^k[ne{T)-n,{T)] 



Where 



(23) 



ne{T) 



1 



ne{T) = 



1 



(24) 



The number density of neutron and proton is defined in Eqns.l2 and 13. The extra 
terms which must be added to energy density and pressure (given in eqns. 14 and 15) 
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Figure 11: Critical temperature versus proton fraction i/p. 



are respectively 



and 
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J^^ dkky{k^ + me^)[ne{T) + ne{T)] 



j_ /-o" dkk\n,{T)+n,{T)\ 



(25) 



(26) 



For a given value of n^, the equations[10,12,13,17,23] arc to be solved self-consistently 
using Eqns.19 and 20 (for fixed S=0,1 or 2) to get /Xp, /Xg, n^, y, T and then can be sub- 
stituted to get pressure and energy density. After getting pressure as function of energy 
density, the TOV equations can be integrated using proper boundary conditions [13] to 
get mass and radius of star at fixed entropy per baryon S. 

Pressure as a function of number density for S = 0, 1 and 2 is shown in Fig. 12. 



One can mark from the figure that the EOS becomes softer as the entropy per baryon 
decreases from S=2 to S=0 which indicates the lowering of mass and radius as shown 
in Fig. 13 and Fig.l4 respectively. For different values of S, the radius, energy density, 
pressure, number density and temperature corresponding to maximum mass is given 
in the following table. 

Table 

Star properties for matter in beta-equilibrium at finite entropy. 



s 


Mmax 
Mq 


R 

(km) 




ric 
(fm-3) 


Pc 

( MeV\ 
\ fmS ) 


T 
(MeV) 





2.18 


12.14 


1230 


0.97 


304.71 


0.0 


1 


2.21 


12.23 


1190 


0.94 


294.85 


27.85 


2 


2.33 


12.45 


1092 


0.85 


272.24 


62.12 



The results in the table reflects the influence of entropy on the gross properties of 
stars. Increase of maximum mass and radius upto S=2, amounts to only a few percent 
of cold star and the temperature of the core is upto 62 MeV which is in fair agreement 
with the results obtained in Ref.[2]. It also shows that for S=2, the thermal contribution 
is larger which results in an increase of mass of about 0.15 Mq compared with that 
for S=0. In neutron stars, the pressure is supported by strongly interacting baryons 
which have smaller contributions to the pressure which causes a smaller increase in 
maximum mass. Thus the compositional variable of EOS has more importance than 
the temperature for the structure of neutron star [2]. But for a White Dwarf which have 
highly ideal EOS [23], the structure and properties are very sensitivity to the entropy 
per baryon. 

The relative concentrations of chemical potentials of n,p and e in beta-equilibrium at 
a fixed entropy per baryon S=l, is shown in Fig. 15. It is clear that /^e increases linearly 
with number density whereas and /ip first decreases and then increases linearly. The 
increase of electron chemical potential with number density implies the abundance of 
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Figure 12: Pressure as a function of number density at fixed entropy per baryon. 

negatively charged particle (electron) which shows that the system is proton rich over 
an extended region of density. It may be seen from the figure that in the very lower 
density region the proton abundance is more, then decreases to some extent and then 
increases linearly in high density region. 

The temperature as a function of energy density at fixed entropy per baryon is 
shown in Fig. 16. The temperature of the star increases for both S=l and S=2 from 
which one can get the critical temperature corresponding to the maximum mass of the 
star. The temperature is maximum at the center of the star (where central energy 
density is about 1100-1200 MeV/fm^ for a maximum mass star) and decreases with 
decreasing energy density which is faster particularly at lower energy densities. This 
implies that the interior of the star maintains a small variation of temperature but falls 



3 



2.5 



1 .5 



0.5 



S = 
S = 1 
S = 2 




500 



1000 



1500 



2000 



(MeV/fm ) 

Figure 13: Star mass (M/Mq) as a function of central energy density at fixed entropy 
per baryon. 



rapidly towards the surface region as density decreases. 



4 Conclusions 

We have studied here asymmetric nuclear matter at finite temperature using the MCH 
model[13]. We have presented how the efi^ective nucleon mass, energy per baryon, 
entropy density, entropy per baryon and pressure behave as a function of density for 
various temperatures. At zero density we find that this model exhibit a phase transition 
at TRi 235 MeV just as obtained in the Walecka model at Tfa 185 MeV. This model 
exhibits the existence of liquid-gas phase transition in asymmetric nuclear matter and 
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Figure 14: Radius versus star mass at fixed entropy per baryon. 

the critical temperature Tc depends sensitively on the proton fraction yp. decreases 
with the decrease of yp as shown in Fig.ll. The value of critical temperature decreases 
from Tc ~ 17.2 to 11.6 MeV for yp — 0.5 to without p-meson and from Tc ~ 17.2 to 
MeV for yp — 0.5 to 0.02 with p-meson. Hence the addition of p-meson seems to be 
very important to the study of properties of asymmetric nuclear matter as it lowers the 
critical temperature. This also shows that even at zero temperature the system remain 
only in gas phase for neutron matter(yp = 0). At fixed temperature and density the 
pressure of the system increases with decrease of proton fraction (shown in Fig. 7) which 
indicates that during the isotherm liquid-gas phase transition pressure can not remain 
constant for asymmetric nuclear matter. We have also studied the EOS and structure 
of PNS with neutrino free charge neutral matter in beta-equilibrium. We find that as 




Figure 15: Chemical potential versus number density in beta- equilibrium at fixed entropy 
per baryon S—1 

the PNS cools from S = 2 to S = 0, the maximum mass and radius exhibit a slow 
decrease. Thus the influence of entropy per baryon or equivalently the temperature, 
on the structure of PNS is not very sensitive. It is also observed that at finite entropy 
per baryon e.g. S=l, the star is proton rich over an extended region of density. The 
temperature varies slowly in the interior of the star but falls rapidly towards the low 
density surface region and the maximum temperature of the core of the star for S=2 
is about 62 MeV. All these results of PNS are in fair agreement with that obtained in 
Ref.[2]. 
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